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ABSTRACT 

Context. The space-borne missions CoRoT and Kepler have provided a large amount of precise photometric data. Among 
the stars observed, red giants show a rich oscillation pattern that allows their precise characterization. Long-duration 
observations allow for investigating the fine structure of this oscillation pattern 

Aims. A common pattern of oscillation frequency was observed in red giant stars, which corresponds to the second-order 
development of the asymptotic theory. This pattern, called the universal red giant oscillation pattern, describes the 
frequencies of stellar acoustic modes. We aim to investigate the deviations observed from this universal pattern, thereby 
characterizing them in terms of the location of the second ionization zone of helium. We also show how this seismic 
signature depends on stellar evolution. 

Methods. We measured the frequencies of radial modes with a maximum likelihood estimator method, then we identified 
a modulation corresponding to the departure from the universal oscillation pattern. 

Results. We identify the modulation component of the radial mode frequency spacings in more than five hundred red 
giants. The variation in the modulation that we observe at different evolutionary states brings new constraints on the 
interior models for these stars. We also derive an updated form of the universal pattern that accounts for the modulation 
and provides highly precise radial frequencies. 

Key words. Stars: oscillations - Stars: interiors - Stars: evolution - Methods: data analysis 


1. Introduction 

The space-borne photometric missions CoRoT 
(|Baelin et al.l 120061) and Kepler (|Borucki et al.l 1 20101) 
have observed many red giants and led to substan¬ 
tial results. It has been shown that oscillations in red 
giants are similar to those seen i n main-sequence solar - 
like stars dDe Ridder et al.1 120091 : iBedding et al.l I2010D . 
Characterization of the oscillation modes in red-giant 
spectra leads to reliabl e estimations of the mass and the 
radius of the stars (e.g.. iKallinger et aI1l2010D . For many of 
them, it is possible to characterize the stellar core, hence 
provide informati o n on the evolutionar y state of the star 
i|Beck et al.l 120111: IBedding et al.l 1201 ll) and on the core 
rotation ( Beck et, al.ll2012l : iMosser et, al.ll2012a|) . 

Pulsating red-giant stars are characterized by two 
distinct resonant cavities, the core and the envelope. 
Gravity(g) waves only propagate in the radiative core. The 
inner turning point of the acoustic(p) modes is located near 
the external edge of the radiative region, and their outer 
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turning point is near the surface of the star. In a red-giant 
star, the transition region between the two cavities is nar¬ 
row, which leads to possibly efficient coupling between p 
and g waves and give s rise to the so-called mixed modes 
(e.g.. Beck et alJl201lt IMosser et ahll2012bl) . These modes 
behave as acoustic modes in the envelope and gravity modes 
in the core. Because of the nature of g-modes, the radial 
modes cannot couple, but dipolar and quadrupolar modes 
can be strongly coupled. 

Radial solar-like oscillations, as observed in red giant 
spectra, have been identified as pressure modes that result 
from acoustic waves stochastically excited by convection 
in the outer layers of the star. The observed pressure 
mode pattern has been depicted in a canonical form , called 
the universal red-giant oscillation pattern (jMosser et al.l 
1201 ll) . This canonical form describes the regularity of the 
pressure mode pattern characterized by two quantities: the 
frequency z/ max of maximum oscillation and the mean fre¬ 
quency difference Av between c onsecu t ive p ress ure modes 
of same angular degree. From IMosser et al.l (2013), we 
can consider the observed pattern to be the translation of 
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t he se cond-order asymptotic pattern described in iTassonll 
(I1980T) at low radial order. 

It has long been predicted that sound-speed discontinu¬ 
ities exist in star interiors a nd t hat they affect the observed 
solar-like oscillations (e.g.. lGoughlll990H . In red giants, the 
main source of discontinuity is the region where helium 
undergoes its second ionization, as shown bv iMiglio et alJ 
(j2010f) . The sharp sound-speed variations produce a mod¬ 
ulation in the observed oscillations frequencies, which is 
called a glitch. In this paper, we aim to investigate the de¬ 
viations from the universal pattern, which are identified as 
glitches. Glitches related to the second helium io nizati on 
zone w ere first measured in the Sun bv iDziembowski et alJ 
(|1991lf . T heir existence in other solar-lik e pulsators was also 
proposed (|Monteiro fe Thompsonlll998D and meas ured for 
the first time in a red giant bv IMiglio et al.l (l2010h . 

Glitches have been also observed in main-seq uence 
stars showing solar-like oscillations (e.g.. iMazumdar et alJ 
|2012|) . The study of glitches in red giants is intended to 
improve the understanding of their internal structure, 
to evaluate the amount of helium present in their en- 
velopes and to charac terize their evolutionary states (e.g., 
IBroomhall et a.1.112014T ). For example, the amplitude of the 
glitch related to the region of second ionization of helium 
is directly correlated to the amount of helium present in 
the envelope. 

We us e stars with an evolutionary state already deter¬ 
mined by iMosser et al.l (|2014f) in order to investigate the 
behaviour of the glitches as a function of the evolution¬ 
ary status. In Section 2, we detail the extraction of the 
radial frequencies of the stars. Section 3 provides a descrip¬ 
tion of the oscillation pattern used as a global reference. In 
Section 4, we characterize the remaining modulation iden¬ 
tified as glitches. In Section 5, we discuss the relationship 
between the measured parameters and c ompare them to 
predictions made by previous models (e.g ., IBroomhall et alJ 
12014 IChristensen-Dalsgaard et ahl 12014 1. Section 6 is de¬ 
voted to conclusions. 

2. Data analysis method 

2.1. Data set 

We used the long-cadence data from Kepler with the max¬ 
imum available length, up to the quarter Q17, correspond¬ 
ing to 1470 days of photometric observation. Original light 
curves were processed and corrected from phenomenolog¬ 
ical effects such as o utliers, jump and drifts according to 
the method of lGarcfa et al.l (|201lh We used the set of 1142 
stars for which Mosser et al.l (12014) deduced the evolution¬ 
ary status from the identification of the mixed-mode pat¬ 
tern. This allowed us to select oscillation spectra with a 
high signal-to-noise ratio as measu r ed by an a ut ocorre la- 
tion amplitude above 50 (Mosser fc Appourchaux j 12009 *1 

2.2. Radial mode fitting method 

As indicated previously, dipolar and quadrupolar modes 
in red giant stars oscillation spectra have a mixed com¬ 
ponent. They have, then, a complex frequency pattern. 
Consequently, we have focused on radial modes only. 


The universal pattern predicts the positions of modes 
using the so-called large separation which is the average fre¬ 
quency separation between modes of the same degree and 
consecutive order. A first estimate of the frequency position 
of radi al mod es is determined with this universal pattern 
( Mosser et al.ll2011 L This guess is refined by automatically 
locating the nearby local maxima in the smoothed spec¬ 
trum within a range of frequencies equivalent to a tenth of 
the large frequency separation. These local peaks are con¬ 
sidered to be reliable modes when their heights are above 
a threshold value corresponding to the rejection of the H 0 
hypothesis with a confidence level of 99.9%. 

We then consider small frequency windows around 
each separate radial mode. We fit a Lorentzian model 
to these modes using the maximum likelihood esti mator 
techniq ue desc ribed in lToutain fe Appourchauxl (Il994l l and 
lAppourchaux (2014j). The structure is fitted as a number 
of Lorentzians plus a backgro und component (Fig. [l]) fol¬ 
lowing the model of lBarban et al.l (|2010T> 

Q 

P{y) = ^-M(/c,zz) +B(v), (1) 

fc=1 

where Q is the number of modes modelled by a Lorentzian 
function M(k, z/), and B(i/) is the background noise in the 
power spectrum modelled as a linear function of u, i.e a + 
bi^k ■ Sever al alt ernative f orms of the background were tried 
(e.g.. iMathur et aPl201lH . and we found that a change in 
the background component has a negligible impact on the 
measured frequencies of the radial modes. Q may be greater 
than 2 when dipole mixed modes are present or when the 
quadrupole modes are split by rotation. 

AI(k, u) is a Lorentzian profile 



where Hk is the height of the Lorentzian profile, Vk is the 
oscillation mode frequency, and Tfc is the mode linewidth 
(FWHM) with Tk = l/(7rr), t being the mode lifetime. 
This method provides realistic uncertainties on the deter¬ 
mination of the mode frequency. The mean value of the un¬ 
certainties we obtained for our set of stars is about 10 nHz. 
We checked that the modelling of the background does not 
modify the fit of the modes frequencies and found that it 
has a negligible impact on the frequency position measured 
for the fitted modes. The position of the radial modes are 
then checked by eye for verification purposes: a non-radial 
mode could be mistakenly identified as a radial mode. 

3. Glitches inference 

The aim of this section is to describe how we measure acous¬ 
tic glitches. This measurement relies on the precise identi¬ 
fication of radial-mode frequency spacing. Since there is no 
unique definition of those spacings, we have to introduce 
different notations. The mean observed large separation is 
denoted (Az/). A weighted mean large separation over a 
broad range of order was used to measure this parameter. 
The observed difference between individual radial mode fre¬ 
quencies is denoted by A i/(n) and Az/ UP (n) is used for the 
reference values provided by the universal pattern. 
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by the frequency difference between two consecutive radial 
modes. The frequency reference for each difference is taken 
as the mid-point as the data values used to derive the differ¬ 
ence. The use of Eq. to derive the frequency separation 
introduces correlations since every radial mode frequency 
is used for estimating two values of A v{n). The frequency 
separation A v{n) is correlated with A v(n ± 2) but uncor¬ 
related with A v(n ±1). This correlation is later taken into 
account in the error budget. 

The results for one star are shown in the top panel of 
Fig. [2] Two notable features can be noticed: an upward 
trend with increasing frequency which we associate with 
the second-order term in the asymptotic expansion, and a 
modulation that we later attribute to glitches. 

3.3. Curvature 


Fig. 1: Power density spectrum (PSD) of the red clump star 
KIC2303367. The peak around 32.1 pHz is a radial mode and 
the peak on the left is the nearby 1 = 2 mode. The Lorentzian 
fits are shown with the solid green line. 


3.1. Determination of the global large separation (Av) 

Since red giant stars pulsate at medium to low radial 
order (the typical radial orders at which red giant stars 
pulsate are about 10 by comparison with around 20 for 
main-sequence stars), second-order terms in the asymptotic 
expansion need to be considered. 

The autocorrelat ion method described in 
iMosser fe Annourchauxl (12009H provides a first ap¬ 
proximation of the global large separation (Av) and of the 
frequency v m ax of maximum oscillation power. As with 
any method, this automated method is affected by the 
so-called realisation noise due to stochastic excitation of 
short-lived waves (Fig. [[]). A way to lo wer th is noise is 
provided by the universal pattern (IMosser et al.l f20Tlh . In 
practice, the correlation of the observed spectrum with a 
synthetic spectrum allows us to substantially reduce this 
noise. It follows that this autom ated meth o d give s precise 
results, even at low frequency (|Hekker et, al.ll2012lh _ 

Th e universal pattern as introduced bv IMosser et al.1 
(I 20 TI for red giants describes the second-order term of 
the asymptotic expansion with a curvature term a: 

i/ up (n,0) = (ii + £+^(n-n max ) 2 ) (Az/), (3) 

where e is the asymptotic offset and n max is the radial or¬ 
der corresponding to the oscillation of maximum amplitude 
defined by n max = z'max/A v — e. Note that n max is not an 
integer. 

3.2. Determination of the frequency differences Av[n) 

We calculate the local frequency separation by computing 
the frequency differences between consecutive radial modes. 
We derive the local A v(n) from the average over adjacent 
modes, 


M \ Di+1,0 ^n—1,0 /,\ 

n) = - 2 -’ ( 4 ) 

where n is the radial order. At the edges of the measured 
radial modes, we cannot use this equation and replace it 


We can deduce the asymptotic dependence of the large sep¬ 
aration with the frequency from the derivative of Eq. m 
with respect to the radial order n: 

A^ UP (n) = (1 + u{n — n max ))(Aiv). (5) 

According to the universal pattern, the variation of 
A u VP (n) is linear and only depends on the curvature term. 
As this asymptotic development is only valid for smooth 
interiors, inner discontinuities will lead to a departure from 
that pattern. In order to identify any modulation due to 
internal discontinuities, we have to remove the curvature 
term. 

We used the relation 

a = 0.015(Ai/)-°- 32 , (6) 

which provides a consistent glitch-free description of the 
curvatu r e term av eraged over a large number of stars 
(IMosser et al.112013 1. In fact, the curvature observed in an 
individual star is affected by the existence of glitches. This 
is illustrated in Fig. [2] (top), where it is evident that the lin¬ 
ear term is steeper when it is derived locally. Using it would 
yield an overestima te of the asymp t otic s econd-order effect 
(IMosser et al.ll2013h . iHekker et al.1 (1201311 showed that the 
curvature is significantly overestimated when only a small 
number of radial modes are used to derive it. We therefore 
use Eq. © as a reliable formulation to specify a. We also 
tested if a change in the a parameter modifies the results 
obtained for the glitch measurement. We found that it is 
important to use a consistent value of a for all stars but 
the precise formulation is not crucial as long as the changes 
are small (see section 4.4). 


4. Measurements 


4.1. Measuring the glitches for red giants 

The methods used for measuring the modulation due 
to the glitches were mostly developed for main-sequence 


_tiy 

stars dMonteiro et al.1119941: iMonteiro fo Thomnsonl 120051 : 


iHoudek fc Clough 200711 . Two discontinuities are mainly 


considered: the base of the convection zone and the helium 
second ionization region. In red giants, the dominant 
discontinuity is t he se cond helium ionization zone, as 
shown by iMiglio et al.l (|2010f l. This work also pointed 
out that there are too few identified modes to take all 
discontinuities into account. Since the discontinuity of the 
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Fig. 3: Apparent curvature a 0 b s as a function of the large 
separation for each red giant stars. Clump stars are indi¬ 
cated by red diamonds and RGB stars by blue triangles. 
Error bars correspond to the lcr uncertainties. The solid 
black line correspond t o the value of the mean curvature a 
found bv lMosser et al.1 (|2013fl 


We fitted an oscillatory component to the resultant fre¬ 
quency variations obtained after removal of the curvature 
term from the measurements: 


r _ A / A \ I ^7 t(v b 'max) 

Og,obs A{Al') COS | Q(Av') ' 


( 8 ) 


Fig. 2: Top: Variation of the large separation as a function of 
frequency for the star KIC1872166. Error bars correspond 
to the l-cr uncertainties. The general trend is due to the 
second-order term (the curvature) and the modulation is 
due to the presence of acoustic glitches. The red dashed line 
is a linear fit over the oscillations to estimate the curvature 
parameter. The blu e dashed line is the value of the mean 
curvature given bv lMosser et al.l ( 2013 1 Bottom: Variation 
of the large separation S n as a function of frequency af¬ 
ter suppressing the curvature term. The dotted-dashed line 
shows a fit obtained for this star using Eq. (|S|). 


second helium ionization is by far the most important 
one in red giants, we can assu me that there is only one 
modulation component (Fig 2. of lMiglio et al.ll20Toh . Most 
of the methods used to characterize the glitches describe 
them in a complex way with a frequency-dependent 
ampl i tude dMonte iro fc T hompsonll2005l : iHoudek fe Goughl 


2007: iMazumdar et al. 1 201 -l|) . Here, due to the low number 

of modes measured and the large uncertainties, we prefered 
to use that simple fit. 

We compute the difference between the observed local 
large separation and the theoretical local large separation 
predicted by the universal pattern: 


<ig, ob s = A v(n) - Ai z UP (n). 


(7) 


where Q is the period of the oscillation expressed in units 
of (Azz), A is the amplitude of the oscillation in units of 
(Azz) and (f> is the phase of the oscillation centered on i/ max . 
We used a y 2 minimisation to fit the three free parame¬ 
ters. The uncertainties were extracted by the inversion of 
the Hessian matrix constructed with a covariance matrix 
(see Appendix 1X1). This method takes the correlation of the 
frequency difference calculated with Eq. Q into account. 
An example of a fit is shown in the bottom panel of Fig. [2] 

In some of the stars we analysed, the fit was unsuccess¬ 
ful. This happens, for example, when the uncertainties on 
the frequencies are too large. We choose to reject such fits 
when the uncertainties on the period parameter were higher 
than the measured value of the period. 

It is important to note that the uncertainties on zz max 
impact on the uncertainties measured for the phase param¬ 
eter. The error on z/ max corresponds typically to a fifth of 
the (Ai/) value. Therefore, following Eq. ©, an error on 
v max will translate into an error on the phase parameter of 
about 27 t/(5 G). With Q ~ 4, this will quadratically add an 
error of 0.3 rad on the phase. It is worth noting, that this 
value is much smaller than the pattern found in Section 4.4 
for this parameter. 

4.2. Period of the modulation 

We measured the dimensionless period Q for 546 stars. The 
variation of the period as a function of the large separation 
is shown in Fig. [4] The dimensionless period Q shows no 
significant trend with (Ai/) over a large range. The clump 
star values show a higher dispersion than the RGB stars. 
We can however note a second branch of RGB stars with 
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Fig. 4: Top: Dimensionless period Q of the modulation mea¬ 
sured as a function of the large separation. Clump stars are 
indicated by red diamonds and RGB stars by blue triangles. 
Error bars correspond to the ler uncertainties. The dashed 
black line indicates the maximum number of radial modes 
observable in a red giant spectrum. Bottom: Histogram of 
the dimensionless period Q. Clump stars are indicated by 
the red dashed line and RGB stars by the blue line. 


much higher period values. As of now, there is no expla¬ 
nation toward this behaviour. The mean glitch periods are 
Q ~ 3.08±0.65 for RGB and Q ~ 3.83T0.88 for clump stars, 
higher than for RGB stars by approximately 30 %. Such val¬ 
ues are similar t o the periods predicted by the models for 
this kind of star dBroomhall et al.ll2014li . 


The values found for the glitch periods are less than 
the typical number of radial modes measured in red 
giant, which is around 6. As the global large separation 
measurement takes all radial modes frequencies into 
account, we infer that glitches are smoothed out when 
the mean large separation is measured globally. As a 
consequence, the glitches have limited influence on a global 
measurement of the large separation. On the contrary, they 
may considerab ly affe ct local measurements, as shown by 
iKallinger et, al.l ( 2012ft . 


Fig. 5: Dimensionless amplitude A of the modulation mea¬ 
sured as a function of the large separation. Clump stars are 
indicated by red diamonds and RGB stars by blue triangles. 
Error bars correspond to the ler uncertainties. 


4.3. Amplitude of the modulation 

The variation of the dimensionless amplitude A as a func¬ 
tion of the large separation is plotted in Fig. [5] It shows that 
the relative amplitude of the modulation A for RGB stars 
increases when the large separation decreases as a conse¬ 
quence of stellar evolution. Clump stars have roughly con¬ 
sistent amplitudes within the uncertainties. There is some 
evidence that the amplitudes for secondary clump stars are 
higher than for RGB stars. We checked that these phenom¬ 
ena were not due to the exclusion of a frequency dependent 
amplitude in the fit. 

We can now derive a relationship between the relative 
amplitude and (Ai/). 


A = a{ Az/) b , (9) 

with a = 0.06 ± 0.01 and b = —0.88 ± 0.05 for RGB stars 
with (Ai/) in /rHz. The absolute amplitude (i.e. the product 
of A and (A v)) shows almost no dependence with (Az/) 
since th e power index b is clo se to —1, as derived from 
models dBroomhall et al.ll2014li . 

Stellar evolution models show that the amplitude of the 
discontinuity in the adiabatic exponent increases when the 
star evolves up the RGB. Such an inc rease is also present in 
the sound-speed profile (e.g.. iHekker et al1l201l b The am¬ 
plitude of the discontinuity strongly depends on the amount 
of helium in the convective envelope of the star. During the 
red giant evolution, a structure change at, e.g., the luminos¬ 
ity bump could bring additional helium in the convective 
envelope of the star. As a result, a larger discontinuity due 
to the helium second ionization zone is expected. It could 
explain why clump stars have a slightly larger observed 
glitch amplitude than RGB stars below the bump. 

4.4. Phase of the modulation 

We now turn to a consideration of the measured phase <f> 
of the oscillation due to the glitch. This parameter shows 
complex variation (Fig. [6]). All clump stars have a phase 
between —2 rad and 1 rad, whereas RGB stars have <j> 
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Fig. 6: Phase <t> of the modulation measured as a function 
of the global large separation. Clump stars are indicated by 
red diamonds and RGB stars by blue triangles. Error bars 
correspond to the la uncertainties. Bottom: Histogram of 
the phase 0. Clump stars are indicated by the red dashed 
line and RGB stars by the blue line. 


of a va riation in internal structure related to the luminosity 
bump dLaearde et al.ll2012i) . 

4.5. On simulated data 

We have performed tests to evaluate the statistical stability 
of the fit of the glitches (Eq. [BJ). First, we have tested the 
influence of pure noise on the universal oscillation pattern. 
Synthetic radial mode frequencies have been created, based 
on the asymptotic pattern and including a sinus-like mod¬ 
ulation. A random shift has been added to each frequency, 
corresponding to a gaussian noise component. We then de¬ 
duced the local large separation from these frequencies. We 
have used the fitting program to extract the components 
of a spurious glitch signal introduced by this noise compo¬ 
nent. The results obtained with one thousand simulations 
indicate that there is a preferance in favour of the detection 
of spurious, low period signals (Q around 2.5) independent 
of the large separation or the signal-to-noise ratio. These 
low periods are associated with random phases and ampli¬ 
tudes. These characteristics are typical of the white noise 
and do not correspond to the signal we observe. Hence, this 
first set of simulation proves that what we observe for the 
phase and period parameters is definitely not due to noise. 

We performed a second test, with the addition of a sinu¬ 
soidal modulation to the oscillation pattern. Various levels 
of noise were considered. The parameters of the sinusoidal 
component could be retrieved, even at low signal-to-noise 
ratio (S/N = 1). These results confirmed that the analysis 
is relevant, even when only five radial orders are considered. 
We however noticed a small bias on the inferred period. 
Further tests have been performed to provide an estimate 
of this bias. For a modulation period G of the sinusoidal 
component, the fit provides a period G* = 1.03G + 0.64. 
Periods were calibrated to take that phenomenon into ac¬ 
count. Similar tests have shown that the phases and am¬ 
plitudes are reliably measured. We also tested for the im¬ 
portance of the number of orders included in the data. We 
found that a modification from 10 to 5 radial orders give 
consistent results within 10 %. 


between 1 rad and 3 rad. This phase shift between clump 
and RGB stars is systematically observed. We investigated 
the influence of the curvature term on these results and 
found that a curvature change as high as 50% does not 
influence the trend. Thus, the difference between clump 
and RGB stars is real and is not an artefact of the data 
processing. 

To investigate the consequences of this phase difference, 
we consider the phase at the order corresponding to the in¬ 
dex n max of the maximum oscillation signal. Following Eq. 
©, if the star has a phase ~ 0 like clump stars, the local 
large separation measured will be overestimated. On the 
contrary, if the star has a phase ~ 7r like RGB stars, the 
local large separation will be slightly underestimated, as 
shown in Section 5.5. This property can be used to deter¬ 
mine the evolutionary stage of the stars, as discussed in the 
section 5.3. 

Additionally, we note the transition to lower values of 
the (f> parameter for RGB stars with a large separation be¬ 
low 6/xHz. This could possibly correspond to the signature 


5. Discussion 

In this section we interpret our results in terms of phys¬ 
ical quantities and compare them to stellar models. The 
modelling of the glitch is then combined with the universal 
pattern to provide a precise fit of radial frequencies. 


5.1. Acoustic radius of the Hell ionization region 

Any structure discontinuity induces a seis mi c sign ature 
characterized by its acoustic radius fe.g.. lGoughlll990H . The 
measured per iod G is directly related to the acoustic depth 
of the glitch (iMazumdar et al.ll2014D 


1 

2G(Au) ~ Tg 



( 10 ) 


where c s is the adiabatic sound speed, r g is the acoustic 
depth of the glitch, r g is its radial position and R s is the 
seismic radius of the star. If the centre of the star is taken 
as a reference, we can define the acoustic radius as: 
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Fig. 7: Acoustic radius of the discontinuity related to the 
second helium ionization zone as a function of the global 
large separation. Clump stars are indicated by red dia¬ 
monds and RGB stars by blue triangles. Error bars corre¬ 
spond to the ler uncertainties. The light green line indicates 
the theoretical acoustic radius of the second helium ioniza¬ 
tion zone for a 1M@ star during the RGB phase. The dark 
green line gives the same information for a 1.4M 0 star. 


T n 



( 11 ) 


We can relate the two physical parameters by the relation 
T g = Tq — T gi with To the total acoustic radius of the star. 
The large separation is also related to the stellar acoustic di¬ 
ameter. However, the global large separation {Av) we mea¬ 
sure is different from the theoretical l arge separa tion, called 
asymptotic large separation Az/ as dMosser et al.Il2013h 



where c s is the sound speed. However, iBelkacem et alJ 
(2013) have shown that {Av) provides a valuable approxi¬ 
mation for deducing the mass up to 10 % and radius of the 
star up to 5 %. Hence, we use it in this work. 

Since the large separation depends on the total acoustic 
radius of the star, we can write T 0 = 1/{2{Av)). Thus, we 
can relate the acoustic radius to the period measured: 


Tg 2{Av) ( X 2 g) ' (13) 

The measured T g are compared to results from CESAM2k 
models in Fig . [71 T he models used are described in 
iBelkacem et al. ( 2012) . The conclusions obtained for these 
m odels are confirm ed by other simulations done by, e.g., 
iMiglio et al.l (l2010il . 


There is a reasonable agreement between the observed 
and modelled values of the acoustic radius for RGB stars 
but a clear discrepancy for clump stars (Fig. [7]). Simulations 
show that a clump star and a RGB star with the same mass 
and radius should have about the same acoustic radius T g . 
This is not what we observe. These are the first extensive 
masurements of the acoustic radius of the helium ionization 


zone in clump stars, the discrepancy that we observe could 
either be due to the difficulty to model clump stars, or to 
another poorly known phenomenon like the mass loss at 
the tip of the red giant branch. We leave an exploration of 
the reasons for this discrepancy to future work. 

5.2. Measuring the helium abundance with glitches? 

iBroomhall et al.l (|2014f ) have shown that red giants with 
similar ( Av) but with different amount of helium show very 
similar glitch amplitude. Owing to the large uncertainties 
we observe on these amplitudes, we have not been able to 
extract information on the helium abundance present in 
these stars. Furthermore, at this stage we do not identify 
any way to improve the measurement and derive quantita¬ 
tive information on the helium content in red giants from 
the glitches. 


5.3. Measuring the evolutionary stage 


Previous work has shown a partition of red giants, depend¬ 
ing o n their evolutionary states, based on the e offset (Fig. 
3b of iBedding et al.l 1201 lb IKallinger et al.l l2012h . This pa¬ 
rameter £ in the asym ptoti c ex pansi on is equal to 1 /4 in the 
asymptotic relation ([Tassoull ll98dH . Its value transferred 
in the universal pattern (Eq. [3j) is significantly different 
and varies with {Av). The link between e as = 1/4 and 
the observed values e((Av)) is however totally explained 
by th e use of differ ent definitions of the large separation 
(IMosser et alJ 120131) . The difference As = £ c ium P — £rgb 
between clump stars and RGB star s with similar (Av) is, 
according to iKallinger et all ([2012T) . —0.15. These results 
are derived from a local analysis of the radial spectrum. 

Our approach, based on a global description with a 
global measurement of the large separation and a generic 
description of the radial oscillation pattern, gives an in¬ 
terpretation of the physical origin of the difference. We 
did not find any e offset between clump stars and RGB 
stars but we observed that the glitch component modi¬ 
fies the local measurement of the large separation with a 
relative variation of —0.5% for RGB stars and +1% for 
clump stars. This translates into a change in the e param¬ 
eter (fe = — (n + e) <51og(Az/) according to the derivation 
of Eq. ©), corresponding to £ c ium P = +0.05 for RGB stars 
and £rgb = —0.1 for clump stars with (Av) = 4/itHz. The 
difference between clump stars and RGB stars is there- 
fore A £ = —0.15 in agreement with local measurements 
(IKallinger et al.ll2012l) . 

The difference reported by IKallinger et al.l (j2012T ) has 
been noted for a vast majority of red giants. We can there¬ 
fore conclude that our results, reduced to a limited subset 
of red giants showing enough oscillation modes, can be ex¬ 
tended to all red giants. Therefo r e, this work justifies the 
analysis made bv IKallinger et al. . (12012) for distinguishing 
RGB and clump stars. Alternatively, measuring the phase 
shift (Eq. 0) provides similar information. 

So, we have now two methods for measuring the 
evolutionary stage of red giants, either based on the m ixed 
modes (e.g.. iBedding et aH 12011 : Mosser et al.l l 2012bh . or 
based on the glitch phases (IKallinger et~ . 20121 . Grosjean 
et al. (2014) recently confirmed that gravity-dominated 
mixed modes in low-mass stars with Av less than 6 ^Hz 
have small heights. This hampers the use of mixed modes 
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Fig. 8: Stellar masses in function of the phase of the modula¬ 
tion. Clump stars are indicated by red diamonds and RGB 
stars by blue triangles. The black box in the right corner 
of the figure corresponds to the typical ler uncertainties for 
these parameters 


for determining the evolutionary stage and reinforces the 
importance of the measurement of the evolutionary states 
only from radial modes. 


5.4. Mass dependence 


Based on the measured (A v) and and on the temper¬ 
ature of the star given by iHuber et al.1 : 201 l ':, we deduced 
the approximate mass and radius of the star from the usua l 
scaling relations dKallinger et al.ll2010l : lMosser et al.ll2013lf . 
We then investigated the mass dependence of the differ¬ 
ent glitch modulation parameters. Only the phase shows a 
clear variation with the stellar mass (Fig. [H]) - We did not 
find any correlations between the mass and other parame¬ 
ters (period and amplitude). Phases of the clump stars have 
a clear mass dependence: clump stars with a high mass have 
a higher phase than clump stars with a lower mass. This 
mass dependence is not seen for RGB stars. 


The observed relation between the phase and the mass 
is clear but remains purely empirical. Its physical ba¬ 
sis needs to be established. More theoretical work is 
theref or e needed, as a follow up to , e.g., iProvost et al.1 
(1199.811: Christensen-Dalsgaard et al.1 (1201411 for p modes 
and Miglio et al.1 ( 2008^ for g modes. 


5.5. Revised asymptotic expansion 

We can now propose a complement to the universal pattern, 
which includes the modulation in order to investigate the 
consequence of the parameter variations in the observed 
shift. To achieve that, we used Eq. ©, corresponding to the 
asymptotic universal pattern for radial modes, and added 
the glitch component modelled in Eq. ©. 

For radial modes, with n max = zz max /(Ai/) — e, we 

get: 


Fig. 9: Synthetic echelle diagram of the radial modes of a 
typical red giant star with a large separation of 3.90 pHz. 
The glitch contributions are overplotted with a dashed red 
line and red diamonds for a typical clump star, and with a 
dot-dashed blue line and blue triangles for a typical RGB 
star. The horizontal dashed black line with a n value of 
approximately 7.6 represents the n max value. Only six ra¬ 
dial inodes are shown in this figure, which is typically the 
number of observed modes. 


Table 1: Asymptotic parameters 


parameter RGB 


clump 


£ 

0.601 +0.632 log(A^) 

^max 

^max/ £ 

a 

0.015(A^)-°' 32 

A 

0.06(A^)-°' 88±O 05 0.07<Az/)-°- 74±0 ' 05 

G 

3.08 + 0.65 3.83 + 0.88 

4> 

1.71 + 0.77 -0.43 + 0.66 


with Af in /j,Hz. 


( OL , •. o 

^ n ~ £ H" ^ max ) 

AG . f 27 r(n - n max ) \\ 

+^sml - - -+ 0ll(A v). (14) 

From lProvost et aTl (Il998h . we derive that the expression is 
in fact purely asymptotic, with the modulation term rep¬ 
resenting the signature of the structure discontinuity. The 
difference compared to the unmodulated pattern is shown 
in Fig.|U As seen above, a local measurement of Av around 
I'max will slightly overestimate the (A v) value for clump 
stars and slightly underestimate it for RGB stars. 

The identification of the modulation provides simultane¬ 
ously a justification and an improveme nt of t he asymptoti c 
red giant universal oscilla tion pattern dMosser et al.ll2f)llll . 
Independent of this work, iMosser et al.l (|2013f ) have shown 
that the e(A^) function observed in the red giant regime is 
the observational counterpart of the second-order asymp¬ 
totic expansion. This work shows that the combination of 
the universal pattern and the modulation (Eq. 1141) provides 
a very precise determination of the eigenfrequency pattern 
with the parameters defined in Table 1. 
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5.6. Curvature and evolutionary status 

Finally, we come back to the influence of the modulation 
on the apparent gradient of the radial frequency separation 
(Eq. [5]). This gradient trans lates into an appar ent curva¬ 
ture in the echelle diagram dMosser et al.ll2013lh Figure [9] 
shows that the curvature of the radial ridge depends on the 
evolutionary status: around ir max , clump stars often show a 
higher curvature than RGB stars. This reproduces the ob¬ 
servations (Fig. [3]) and justifies the use of a mean curvature 
term (Eq. 0 used for removing the second-order effect. It 
demonstrates the relevance of the global picture provided 
by the universal pattern with glitches (Eq. [T-IT) for describ¬ 
ing the red giant radial oscillation pattern. 


6. Conclusion 

In this study, we measured the frequency position of radial 
modes very precisely using a peakbagging method. From 
these frequencies, we calculated a mean value of A v. Using 
the universal red giant oscillation pattern, we analysed the 
variation of A v as a function of the frequency. Modulations 
could be measured for half of the data set. We character¬ 
ized the modulation in more than five hundred red giants 
as a glitch signature due to the second helium ionization 
zone. A value of the acoustic depth of this discontinuity 
has been extracted from this modulation and found consis¬ 
tent with predictions for RGB stars but not for clump stars. 
A different behaviour for the modulation was found in the 
phase parameter between the different evolution stages of 
the stars which leads to an interpretation of the results 
found bv iKallinger et all (l2012h . The mass dependence of 
the phase of the modulation in clump stars needs some 
more theoretical work to be fully understood. Even if ex¬ 
tracting any direct information on the helium content is 
certainly very difficult, this work opens the way for testing 
the second helium ionization zone in a large set of stars with 
ensemble asteroseismology. Finally, it shows that the global 
measurement of the large separation is less biases than the 
local measurement, hence provides better estimates of the 
stellar mass and radius. 


Appendix A: Estimation of the uncertainties 

The x 2 method is used here to fit a model (j/j(a), i = 1 ,N) 
with a set of parameters a over a range of data points yi 
(which correspond to a combination of frequencies). It is 
written 


X 


2 



(A.l) 


with Oi the error on each data point. 

The uncertainties on the parameters can be retrieved 
by the inversion of the Hessian matrix 


Hu = 


d\ 2 


N 


da k dai 




1 dyj( a) dyj( a) 
of da k dai 


(A.2) 


when data are not correlated. The second-derivative terms 
were neglected. 


Since the frequencies considered are correlated, the pre¬ 
vious relations can not be valid. Another definition of Eq. 
ED must be used, which includes a covariance matrix tak¬ 
ing the correlation between the parameters into account. 
Second derivative terms are usually neglected. 

The covariance matrix indicates the correlation between 
the different data points. The data points are calculated 
using the set of radial mode frequencies. If we take two 
combinations of radial frequencies (/ and g ): 

/ = E u n v n (A.3) 

n 

g = '^2,v n v n (A.4) 

n 

where v n are the radial frequencies, u n and v n are the co¬ 
efficients used to calculate the large separation with the set 
of frequencies. For example, if / correspond to the calcula¬ 
tion of the local large separation A v{n) = {v n +1 — v n - 1 )/ 2 , 
we then have u n -\ = —1/2, u n +1 = +1/2 and the other 
coefficient u n equal to zero. 

The elements of the covariance matrix C are then cal¬ 
culated as 


Cij — U. rl: Vri Oy .j , /X . 5) 

n 

The elements of the covariance matrix are then set to zero 
when there is no correlation between the data points. 

The x 2 calculation is then modified following 

X 2 = (y-y(a)) T W(y-y(a)), (A.6) 

with W = C” 1 the inverse of the covariance matrix. If 
there is no correlation between the data points, we retrieve 
Eq. ED. The Hessian matrix calculation has to take the 
correlation between the parameters into account as well: 


r dx 2 
oa k oai 

f , r f dyjja) d yi (a) , d yi {&) d yj {a) \ 

ee dcLk dai + dak da ) 

9—1 ■’1 = 1 ' ' 


(A.7) 


Second-derivative terms were neglected like in Eq. El 
The errors are then extracted from the eigenvalues of 
the inverse of this Hessian matrix. 
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